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Abstract 

We consider a leptophobic Z' scenario in a flipped SU(5) grand unified theory obtained from heterotic 
string theory. We show that the allowed Z' mass, flavor conserving and flavor changing couplings of 
the Z' to the down-type quarks are strongly constrained by the mass difference in B s — B s system 
and the four branching ratios of B — > ttK decays. It is shown that even under these constraints large 
deviations in direct and/or indirect CP asymmetries of B — > ttK decays from the SM expectations 
are allowed. Especially it is possible to accommodate the apparent puzzling data in B — » ttK CP 
asymmetries. 
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I. INTRODUCTION 



The flavor changing neutral current (FCNC), in particular the b — > s transition, is a sensitive 
probe of new physics (NP) beyond the standard model (SM) of particle physics. The processes, 
such as B -> X s7 [1], B -> tt/T [2], 5 -> [3], 5 -> [4], and B s -> ^+^~ [5] which 

are dominated by the 6 — > s transiton have attracted much interest because they still allow 
much room for large NP contributions. The experimental data for some of them show apparent 
deviations from the SM predictions [2-4]. 

A viable NP scenario which may give large contribution to b — > s FCNC is leptophobic 
Z' [6]. Extra U(l) gauge groups appear naturally in many extensions of the SM. If some of 
them remain unbroken down to the electroweak scale, the Z' can be light and affect low energy 
phenomenology. In addition, there is a possibility that the new neutral gauge boson does not 
couple to leptons. This kind of Z' gauge boson is called leptophobic. In case the Z' does not 
mix with the SM Z boson, the strong constraints from the electroweak precision tests can be 
avoided. Explicit leptophobic Z' model with these properties has been constructed by Lopez, 
Nanopoulos and Yuan [7] in heterotic string theory. 

The model has gauge group, G = G ohs x Chidden x <-nj(i), where G ob s = SU(5) x U(l), 
Chidden = SU(4) x SO(10), and Gu(i) = U(l) 5 . It also has 63 massless matter fields. It can 
be shown that the Z' gauge boson can be light and leptophobic without mixing with Z. The 
additional feature of this model is that the Z' coupling is generation dependent. Therefore 
tree-level FCNC is generated in general. Since u c and L belong to the same multiplet which do 
not couple to the leptophobic Z 1 ', the Z' coupling to -u-quarks maximally violates parity. And 
the disparity of the couplings to the d c and u c provides additional source of isospin breaking. 

We study the contribution of leptophobic Z' to B® — B° s mixing and non-leptonic decays 
B — > ^K's with this model in mind. However, our analysis can be easily extended to other 
(leptophobic) Z' models allowing tree-level FCNC's. 

The measurement of the mass difference, Am s , in the B° s —B° s system CDF [8] collaborations 

Am? p = 17.77 ± O.lO(stat) ± 0.07(syst) ps" 1 (1) 
is consistent with the SM calculations [9] 

Amf M = 22.57+H ps" 1 . (2) 
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This constrains many NP models [10-13] including Z 1 models, MSSM models, etc. In [6], we 
showed that the Am s constraint on the leptophobic Z' is much stronger than the previously 
considered one [14] from the semi-leptonic -B-decays. In this paper we extend the analysis in [6] 



to include the case where Z' couples to both left-handed and right-handed quarks simultane- 
ously. When both couplings exist at the same time, we will see that the Am s constraint is not 
enough to set the upper bound on the sizes of the FCNC couplings. We go beyond the [6] and 
demonstrate that the additional constraints are available, i.e. the four branching ratios (BR) 
of B — > irK's, which can give the limits even in the simultaneous existence of both left- and 
right-handed couplings. 

The charmless non-leptonic decays B — > irK have been measured precisely enough to probe 
the electroweak amplitudes [2]. The experimental data indicate that while the BRs are con- 
sistent with the SM expectations, some of direct and indirect CP asymmetries show apparent 
(but still debatable) deviations from the SM [15]. Accepting this discrepancy seriously, we 
can see the electroweak penguin sector is the best place to search for NP [15]. We will see 
that in our model the leptophobic Z' can give large contributions to the electroweak penguin 
amplitudes while satisfying the Am s and the BR(B — > nKys. In addition we will see that the 
predicted direct and indirect CP asymmetries can accommodate the discrepancies between the 
SM predictions and measurements simultaneously. 

We note that the merit of our model in explaining the data comes from (i) it automatically 
evades the stringent constraints involving leptons, such as LEP I data, B s — > (ii) there 

are new CP violating phases, and (iii) the characteristic isospin breaking interaction in this 
model can generate large electroweak penguins. 

The paper is organized as follows: In section II we briefly describe our model, the leptophobic 
Z' model in the stringy flipped SU(5) theory. In section III we calculate the Z' contribution 
to the Am s when Z' couplings to the quarks have both handedness. In section IV we consider 
the constraints imposed by the BR{B — > irK) and predict the deviations in the direct and 
indirect CP asymmetries from the SM expectation and compare with the experimental results. 
We conclude in section V. 

II. THE MODEL 

The leptophobic Z' model can occur naturally in some grand unified theories (GUTs) and 
string theories. It naturally avoids stringent low energy constraints thanks to the absence of 
couplings to charged leptons and light neutrinos. There are at least two known mechanisms 
that can generate the leptophobic Z 1 . The first one is obtained via dynamical mixing between 
the U(l) and U(l)' in the E% GUT [16]. The other scenario of leptophobia is obtained in the 
stringy flipped SU(5) GUT in the heterotic string theory [7]. In this paper we consider only 



the latter scenario only because it is more relevant to the B — > 7r.fr decays. 

In the flipped SU(5), the SM particles appear in three copies of the representations 



F = (10, l -) = {Q, d c , z/}, f = (5, - |) = {L, u c }, 



X 



{e c }. 



(3) 



The new neutral gauge boson Z' can be leptophobic if it does not couple to 5 and 1, while the 
quarks in 10 still couple to it [7]. 

In addition to its own beauty this scenario has the following phenomenologically interesting 
features: 

• The new Z' coupling is generation dependent and can generate FCNC processes. 

• The FCNC couplings allow large CP violation. 

• It violates the isospin symmetry in the right-handed up- and down-quarks. 

• The new gauge boson interaction maximally violates the parity in the up-quark sector. 
In the mass eigenstates the interactions of Z' gauge boson with the quarks can be written as 
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u + d^P L 



d + dYPR 



d , 



(4) 



where 5 parameterizes the size of the new gauge coupling relative to the SM coupling and 
is expected to be of 0(1). The c = diag(ci, C2, C3) represent the generation-dependent U(l)' 
quantum numbers [7], and Vl <R (q = u,d) are unitary matrices diagonalizing the quark mass 
matrices. The explicit sets of values for Cj(i = 1,2,3) derived from heterotic string theory are 
given in [7]. Since V*, Vr are unknown, we do not take specific values of q's in [7] and take 
them as free parameters. 

We introduce complex parameters, L and R, 
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to represent the b — > s FCNC couplings. For comparison with [6], we kept the factor 2 in (5). 
Then 
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L%s LllJ b L Z'» + R z sb s Rlll b R Z'» 



h.c, 



(6) 
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where 5 in (4) is absorbed into L and R. To calculate the B — > ttK decay amplitudes we also 
need the Z' couplings to the first generation quarks. Although Q and d c have the same U'(l) 
charges, the mixing effect in (4) can give different couplings to the Z' in general 



C(Z'qq) = 



92_ 

cos 9 



w 



-5 Z'» [u lfl c u L P L u + d lfl (c d L P L + c d R P R )d] , 



(7) 



where we defined 
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j 11 
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j 11 
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From the structure of CKM matrix, we assume the couplings to the left-handed quarks are 



approximately equal, i.e. c\ 



c|. However, in general c d R can be different from c\. 



As mentioned above, the absence of c R is the characteristic feature of the leptophobic flipped 
SU(5) scenario. Note that since 5 and c 9 L are unknown, c 9 L can always be absorbed to 5. In 
addition, 5 does not appear in the expression for Am s and it can absorbed into Lf fe ' or B^ h in 
the B — > tiK amplitudes. So we fix 5 = c 9 L = 1 from now on. 



III. B° - B s MIXING 

In general the Z' can couple to both left- and right-handed quarks simultaneously as can 
be seen in (4). Then we need to extend the operator basis beyond the SM one in the effective 
Hamiltonian describing B® — B° s mixing. 

The most general AB = AS = 2 process is described by the effective Hamiltonian [17]: 

5 3 
i=l i=l 

where 

Qi = n-i^Af^L 

Q 2 = s R bls R bi 

q 3 = s a X4bi 

Q 4 = s a R b a L f L b R 

Qs = s a A4b R (10) 

and the operators <5i,2,3 are obtained from the Qi,2,3 by the exchange L <-> R. Here q R ^ = 
Pr,l Q, with P R L — (1 ± 7s)/2, and a and (5 are color indices. 

In our model the nonvanishing Wilson coefficients at Mz> scale are simply given by 



Here Q5 is the additionally generated operator compared with [6]. The renormalization group 
running down to mj scale mixes Q$ with Q4 and we get [17] 



dGO ~ 0.801 C X (M Z ,), 
CM ~ 0.697 C 5 (M Z ,), 

cm ~ 0.886 c 5 (m z 0- 



The other operators are not generated at all and we get C 2 (n b ) = C 3 (/jL b ) = 0. 



Now we can calculate the B® — B s mixing matrix element 

M s l2 = M*i M + M[f = M$ M (l + R), 



where R = M'f/M'f M . The SM contribution M^ M is given by [18], 

^2 JI/T2 



-,SM 



M s,s m = ^m Bs (/ Bs <)VoK) (W/J 2 , 



where 5 Bs ~ Bi^fa,^)] -6 / 23 [1 + 1.627a a (^ 6 )/(47r)]. The Z' contribution 
M[f = l -M Bs f Ba [(CM + CM)BM 
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+ - 



M 
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involves additional hadronic parameters, B^(fi b ) and B 5 (/i b ). 

The mass difference in the B° s —lf s system, Am s is obtained by 



Am s = 2\M{ 2 \ 



In the SM, we get 



Amf = (22.5 ±5.5) ps" 1 , 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



where the nonperturbative hadronic parameters fs s and Bb s are the main sources of the un- 
certainty. We used the value 

f Bs B l £ = (0.295 ±0.036) GeV, (18) 

(HP+JL)QCD 

which is the combined lattice result [9] from JLQCD and HPQCD. For other parameters, we 
used a s (fi b ) = 0.22, rj B = 0.551, mf^(m t ) = 162.3 GeV and V ts = 0.04113 [19]. 

In Figure 1 the allowed region by Am s alone is shown. We fixed mz> = 700GeV which is 
above the experimental lower bound [20] and scanned the weak phases (/>f^ = arg(Lf b '(i?f b ')) 
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FIG. 1: The allowed region in ( 
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plane by Am s alone. 



from to 27r independently. We can see that the sizes of FCNC couplings, \L%\ and \R%\, arc 
restricted typically to be less than ~ 0.1 for most values of the scanned parameters, which is 
consistent with [6] . However, the two additional bands appear in this case due to the cancelation 
between Lf fe ' and R%,- These regions extend indefinitely and show that the Am s alone is not 
enough to constrain both and R% simultaneously. We will show that the BR(B — > irK) 
can give upper bounds on both |L^,'|' S and \R%\ in the next section. 



IV. B irK DECAYS 

The B — > 7i K decays are dominated by the b — > s QCD penguin diagrams. The subdomi- 
nant electroweak penguin contribution is also sizable and may play important role in probing 
the NP as mentioned in the Introduction. The current experimental data in Table I show the 
branching ratios are quite precisely measured and the so-called R c /R n puzzle [21] has disap- 
peared. Therefore we take the four BRs as the additional constraints to the Am s constraint 
considered in section III. As we will see in a moment, they are orthogonal to and as strong as 
Am, constraint. 



Mode BR[10- 6 } A CP S C p 

B+ tt+K 23.1 ± 1.0 0.009 ±0.025 
B + tt°K + 12.9 ± 0.6 0.050 ± 0.025 
B° tt-K + 19.4 ± 0.6 -0.097 ± 0.012 
B°^it K 9.9 ±0.6 -0.14 ±0.11 0.38 ±0.19 



TABLE I: Branching ratios, direct CP asymmetries Aqp, and mixing-induced CP asymmetry Sep (if 
applicable) for the four B — > irK decay modes. The data are taken from Refs. [22] and [23]. 



In the SM the P — > nK decay amplitudes can be written in terms of topological amplitudes: 

A{B+ -> n+K°) = -Pi - l -P[g + 
V2A(B+ - ^K+) = Pi - P EW - 2 -P E ^ - [T' + C + P; c ) 

A(B° - n-K+) = Pi- \p'^ - (V + Pi) e* 
V2A(B° - Tr ^ ) = -p; - P EW - l -P^ - (C - P; c ) e * (19) 

where other small annihilation and exchange amplitudes are neglected. Here the weak phase, 
7, dependence has been explicitly written. The primes denote the b — > s transition. The P t ' c 
(P„ c ) is the QCD penguin amplitude with t, c (u, c) quarks running inside the loop. The tree 
(color-suppressed tree) diagrams are represented by T' (C). The P EW is the electroweak (color- 
suppressed electroweak) penguins and related to the T'(C') by flavor SU(3) symmetry [24]: 

where Cj (i=l,2,9,10) are the Wilson coefficients and R = \V ts V t l/V us V* b \. 

In the SM, from the loop-, color-factor and the hierarchy of CKM matrix elements, we expect 
the following hierarchies: 

o(i) \p; c \, 

o(A) mip EW i, (21) 

0(A 2 ) |C"|, \Puc\, |P E wl> 
0(A 3 ) 

However the experimental data in Table I are not fully consistent with these hierarchies. 
Specifically A CP (B + -> n°K + ) ^ A CP (B° -> 7^+) and ^(P -> vr ^ ) ^ sin 2/3 require 
[C'/T'l = 1.6 ± 0.3 [15]. This large ratio is inconsistent with the SM expectation (21) which is 
supported by theoretical calculations [25-27] . And the implications of this apparent discrepancy 
have been considered in many NP models [2] . This puzzle can be solved most naturally if NP 
is introduced in the electroweak penguin amplitude [15]. 

In this paper, as mentioned above, we use the four PP(P — > 7r.fr) 's to constrain the sizes of 
Z' couplings. Using the remaining parameters after imposing Am s and the PP(P — > 7r.fr) 's, we 
predict the direct and indirect CP asymmetries which show apparent deviations from the SM. 
We show that the predictions for A CP (B + -> tt°K+), A cp (B° -> tt^ + ) and S CP (B° -> 7r°K°) 
in our model are in the right directions to the experimental results. 



To calculate the SM predictions for the BR's and CP asymmetries, we use the NLO cal- 
culations in the perturbative QCD (PQCD) results [26]. In our model the Z' contributions 
can change the amplitudes, -P/ c(nc ), Pew ■ The NP contributions are calculated using the naive 
factorization method and their strong phases are assumed to be equal to the corresponding SM 
diagrams. We will also discuss the effect of the NP strong phase later. The Z' contribution 
to the topological amplitudes are written in terms of the Wilson coefficients in the standard 
operator basis [28] as 

,K„ \ <~ , K- 



P'(Z') = -A t (a 4 + rf a 6 ) - (a 4 + rf a 6 ) 
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Pew( z ) = ij^t 



Pew(Z') = (aio + rf a 8 ) - (a 10 + rf a 8 ) A wK e ld ™, (22) 
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where X t = V t * s V t b, a, = Cj + Cj±i/3 (+(-) for odd (even) i), rf = 2m K /m b (m s + m 9 ) (with 
m g = (m u + m d )/2), A nK ( Kn ) = G F (m 2 B — m1)F^ K ^ f K ^/ y/2, and the <5's are the corresponding 
strong phases obtained in [26]. The Oj's are Wilson coefficients for the chirality flipped operators. 
In our model the Wilson coefficients at M z > scale are 



m% T z>c u L + 2c d L 

■z> 

2 u 

-X t C 5 (M z ,) = 5^fL%^--^, C 6 (M Z ,) = 0, 



\ t C 3 (M z ,) = S^Li L ' \ C 4 (M Z ,) = 0, 
m z , 3 



m z , 

~\C 7 {M Z ,) = 6™fL%(c R - 4), C 8 (M Z ,) = 0, 

Z, 777 g > 

-h t C 9 (M Zf ) = 5™fL*(c u L - c d L ), C W {M Z .) = 0. (23) 

Zi 777 

There are also the chirality flipped operators to the SM operators. Their Wilson coefficients, 
Q's are obtained by exchanging <-> R%, c q L <-> c 9 R . 

In Figure 2, we show the allowed topological amplitudes in (P/ c , P EW ) plane (a) and in 
{P'ewi ^ew) pl ane (b)- We can see that P/ c is strongly constrained to lie in the region (49,52) 
eV by the BR(B — > 7ri^)'s, whereas sizable deviation from the SM predictions are possible for 
Pew an d -Pew- 

As mentioned in the previous section, the Figure 3 shows that, given c q L and c R , the flavor 
changing couplings, |Lf fe '| and \R%\, are constrained by the four BR{B — > 7rfT)'s in addition 
to the Am s . Since the experimental measurements are quite precise now and the theoretical 
calculations have still large errors, we allowed 3-a range for the BRs. For the plot, we set 
M Z ' = 700 GeV and c R = 1,0.5,0 from the left, respectively. Since the BR(B -> ixK) decays 
are most sensitive to the P EW which is maximized at c d R = 1 and vanishes at c d R = 0, the 
constraint is strongest for c d R — 1. 
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FIG. 2: The correlations between P' tc and -Pg W (a) and between -P^^ and i^W C 3 ) f° r = ^00 
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FIG. 3: The allowed region in ( , Rf. ) plane by Am s and the four BR(B — > 7r?T)'s. We fixed 



4 = 1-0, 0.5, 0.0 from the left. 

Now we predict the direct and indirect CP asymmetries using the parameter set allowed by 
Am s and the four BR(B — > 7rK)'s. We are especially interested in the correlation between 
the two direct CP asymmetries, Aqp(B + — > tt°K + ) and Acp(B° — > 7T~K + ), and the indirect 
CP asymmetry, Scp{B° — > tt°K ) because they show the apparent deviations from the SM 
predictions. 

The predictions for the Aqp(B + — > tt°K + ) and Aqp(B° — > tt~K + ) are shown in Figure 4 
for Mz> = 700 GeV. In these figures the errors for the SM predictions are also obtained from 
[26] . For the NP predictions we fixed the SM to the central values and we did not include the 
hadronic uncertainties. Although the SM results are consistent with the experimental data at 
2-cr level, the Z' contribution can accommodate the current data at l-a level for c d R = 1.0, 0.5. 
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FIG. 4: The predictions for .4 C p(£ + -> vr°K+) and A C p(£° -> Tr^iT^) for M Z / = 700 GeV and (a) 
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FIG. 5: The correlation between A C p(B+ -> vr ^) and Scp(B° 
(a) 4 = 1-0. (b) 4 = °- 5 > ( c ) 4 = 



vr ^ ) for Af Z ' = 700 GeV and 



The value c R = 0.0 cannot explain the data. These results are consistent with [15] which claims 
that the current data require large NP contributions at the electroweak penguin sector. 

The predictions for the correlation between Aqp(B + — > tt°K + ) and Scp{B° — > 7i°K°) are 
shown in Figure 5. While it is difficult to get Scp(B° — ► 7r 



as low as ~ 0.38 which is the 



central value for the current experiments, it is possible to accommodate both Scp{B° — > 7i°K°) 
and Acp(B + — > n°K + ) simultaneously for c d R = 1.0,0.5 (Figure 5(a), (b)). Again it is difficult 
to get large deviations from the SM prediction for the Scp(-B — > tt K ) for c d R = 0.0. Therefore 
the value c d R = 0.0, corresponding to P^ W (Z') = 0, is disfavored even if we include the hadronic 
uncertainties in the calculation. 

Until now we fixed the NP strong phases to be equal to the corresponding SM strong phases. 
In general, they may not be equal to each other. To see the effect of the NP strong phases, 
now we allow the strong phase of the dominant NP electroweak penguin 5' EW to take arbitrary 




FIG. 6: The \P^ W \ (a), A C p(B + -► ir°K + ) (b), and S C p{B° -> (c) as a function of strong 

phase, 5^^, of the electroweak penguin. We fixed c R = l,M z > = 700 GeV. 
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FIG. 7: A scattered plot in (M^/,|i?^'|) plane. For this plot we imposed A C p(B + -» 71-°-^+), 
ylcp( J B° -> tt~K + ), and Scp(-B° -> vr°if°) constraints as well as the Am s and BR{B -> vriTl's. 

values. Figure 6 shows that it can give strong impact on the Acp{B + — > tt°K + ). (For these 
plots we fixed c d R = 1, Mz> = 700 GeV.) The reason is that if 5 EW is equal to the strong phase 
of the dominant QCD penguin, the NP electroweak penguin contribution to Aqp{B + — > 7f°K + ) 
vanishes at the leading order of P' EW {Z') / P/ c , independent of the weak phases 4>^ R y The effects 
of 5' EW on l-P^jyl and Sqp(B — > tt°K°) are rather minor, and any value of 5' EW can successfully 
explain the Scp{B° — > tt°K ) anomaly. 

Now we consider the mass dependence of the Z' gauge boson. In Figure 7 we can see that 
the Mz> as large as 5 TeV which is beyond the LHC reach can accommodate the data with 
\R% \ ^ 0.3. (We fixed c d R = 1 for the plot.) The parabolic shape can be understood from (23) 
because the parameter can be fixed in terms of \R%\ by Am s . 



V. CONCLUSIONS 



We considered the leptophobic Z 1 model in the flipped SU(5) GUT obtained from heterotic 
string theory [7]. This is phenomenologically interesting because it contains ingredients which 
can possibly explain the apparent deviations from the SM predictions in the B — > 7r.fr decays: 

• The new Z' coupling is generation dependent and can generate FCNC processes. 

• The FCNC couplings allow large CP violation. 

• The couplings also violate the isospin symmetry and can give large contributions to the 
electroweak penguins, P EW and -Pew- 

We found that if we include the left- and right-handed FCNC couplings Lf fe ' and R% simul- 
taneously, we cannot obtain the absolute upper bounds for them contrary to [6] where it was 
assumed that only a single coupling exists at a time. If we impose the additional constraints, 
BR{B — > 7rfT)'s, with some reasonable assumptions, we can constrain |Lf 6 '| and \R%\- 

We predicted the CP asymmetries, A C r(B + -> n°K + ), A CP (B° -> n~K + ), and S C p(B° -> 
n°K°). Interestingly enough all of them are consistent with the current experimental results 
when the isospin breaking coupling, c R , is non-vanishing. The case for c d R = where there is 
no Z' contribution to the electroweak penguin is disfavored from the current data. 

Acknowledgement The work of C.S.K. was supported in part by CHEP-SRC Program, and in 
part by the Korea Research Foundation Grant funded by the Korean Government (MOEHRD) 
No. KRF-2005-070-C00030. The work was supported in part by the Korea Research Foundation 
Grant funded by the Korean Government (MOEHRD) No. KRF-2007-359-C00009 (SB). J.H.J 
was supported by the National Graduate Science & Technology Scholarship funded by the 
Korean Government. 



[1] A. L. Kagan and M. Neubert, Phys. Rev. D 58, 094012 (1998) [arXiv:hep-ph/9803368]; A. L. Ka- 
gan and M. Neubert, Eur. Phys. J. C 7, 5 (1999) [arXiv:hep-ph/9805303]; S. Baek and P. Ko, 
Phys. Rev. Lett. 83, 488 (1999) [arXiv:hep-ph/9812229]; T. Goto, Y. Y. Keum, T. Nihei, Y. Okada 
and Y. Shimizu, Phys. Lett. B 460, 333 (1999) [arXiv:hep-ph/9812369]. 

[2] R. Fleischer and J. Matias, Phys. Rev. D 61, 074004 (2000) [arXiv:hep-ph/9906274]; J. Ma- 
tias, Phys. Lett. B 520, 131 (2001) [arXiv:hep-ph/0105103]; R. Fleischer and J. Matias, Phys. 
Rev. D 66, 054009 (2002) [arXiv:hep-ph/0204101]; T. Yoshikawa, Phys. Rev. D 68, 054023 



(2003) [arXiv:hep-ph/0306147]; A. J. Buras, R. Fleischer, S. Recksiegel and F. Schwab, Nucl. 
Phys. B 697, 133 (2004) [arXiv:hep-ph/0402112]; V. Barger, C. W. Chiang, P. Langacker and 
H. S. Lee, Phys. Lett. B 598, 218 (2004) [arXiv:hep-ph/0406126]. S. Nandi and A. Kundu, 
arXiv:hep-ph/0407061; S. Mishimaand T. Yoshikawa, Phys. Rev. D 70, 094024 (2004) [arXiv:hep- 
ph/0408090]; Y. Y. Charngand H. n. Li, Phys. Rev. D 71, 014036 (2005) [arXiv:hep-ph/0410005]; 
X. G. He and B. H. J. McKellar, arXiv:hep-ph/0410098; S. Baek, P. Hamel, D. London, A. Datta 
and D. A. Suprun, Phys. Rev. D 71, 057502 (2005) [arXiv:hep-ph/0412086]; Y. L. Wu and 
Y. F. Zhou, Phys. Rev. D 72, 034037 (2005) [arXiv:hep-ph/0503077]; M. Gronau and J. L. Ros- 
ner, Phys. Rev. D 71, 074019 (2005) [arXiv:hep-ph/0503131]; C. S. Kim, S. Oh and C. Yu, Phys. 
Rev. D 72, 074005 (2005) [arXiv:hep-ph/0505060]; S. Khalil, Phys. Rev. D 72, 035007 (2005) 
[arXiv:hep-ph/0505151]; H. n. Li, S. Mishima and A. I. Sanda, Phys. Rev. D 72, 114005 (2005) 
[arXiv:hep-ph/0508041]; R. Arnowitt, B. Dutta, B. Hu and S. Oh, Phys. Lett. B 633, 748 (2006) 
[arXiv:hep-ph/0509233]; Y. D. Yang, R. Wang and G. R. Lu, Phys. Rev. D 73, 015003 (2006) 
[arXiv:hep-ph/0509273]; C. W. Bauer, I. Z. Rothstein and I. W. Stewart, Phys. Rev. D 74, 034010 
(2006) [arXiv:hep-ph/05 10241]; W. S. Hou, M. Nagashima, G. Raz and A. Soddu, JHEP 0609, 
012 (2006) [arXiv:hep-ph/0603097]; S. Baek, JHEP 0607, 025 (2006) [arXiv:hep-ph/0605094]; 
S. Baek and D. London, Phys. Lett. B 653, 249 (2007) [arXiv:hep-ph/0701181]; C. H. Chen, 
C. S. Kim and Y. W. Yoon, arXiv:0801.0895 [hep-ph]; M. Imbeault, S. Baek and D. London, 
arXiv:0802.1175 [hep-ph]. 

C. Dariescu, M. A. Dariescu, N. G. Deshpande and D. K. Ghosh, Phys. Rev. D 69, 112003 (2004) 
[arXiv:hep-ph/0308305]; E. Alvarez, L. N. Epele, D. G. Dumm and A. Szynkman, Phys. Rev. 
D 70, 115014 (2004) [arXiv:hep-ph/0410096]; Y. D. Yang, R. M. Wang and G. R. Lu, Phys. 
Rev. D 72, 015009 (2005) [arXiv:hep-ph/0411211]; P. K. Das and K. C. Yang, Phys. Rev. D 
71, 094002 (2005) [arXiv:hep-ph/0412313]; C. S. Kim and Y. D. Yang, arXiv:hep-ph/0412364; 
S. Baek, A. Datta, P. Hamel, O. F. Hernandez and D. London, Phys. Rev. D 72, 094008 (2005) 
[arXiv:hep-ph/0508149]; C. S. Huang, P. Ko, X. H. Wu and Y. D. Yang, Phys. Rev. D 73, 034026 
(2006) [arXiv:hep-ph/0511129]. 

B. Dutta, C. S. Kim and S. Oh, Phys. Rev. Lett. 90, 011801 (2003) [arXiv:hep-ph/0208226]; 
G. L. Kane, P. Ko, H. b. Wang, C. Kolda, J. h. Park and L. T. Wang, Phys. Rev. D 70, 035015 

(2004) [arXiv:hep-ph/02 12092]; S. Baek, Phys. Rev. D 67, 096004 (2003) [arXiv:hep-ph/0301269]; 
B. Dutta, C. S. Kim, S. Oh and G. h. Zhu, Phys. Lett. B 601, 144 (2004) [arXiv:hep-ph/0312389]; 
M. Endo, M. Kakizaki and M. Yamaguchi, Phys. Lett. B 594, 205 (2004) [arXiv:hep-ph/0403260]; 
J. F. Cheng, C. S. Huang and X. H. Wu, Nucl. Phys. B 701, 54 (2004) [arXiv:hep-ph/0404055]; 



S. Khalil, Mod. Phys. Lett. A 19, 2745 (2004) [Afr. J. Math. Phys. 1, 101 (2004)] [arXiv:hep- 
ph/0411151]. 

[5] C. S. Huang and Q. S. Yan, Phys. Lett. B 442, 209 (1998) [arXiv:hep-ph/9803366]; K. S. Babu 
and C. F. Kolda, Phys. Rev. Lett. 84, 228 (2000) [arXiv:hep-ph/9909476]; P. H. Chankowski 
and L. Slawianowska, Phys. Rev. D 63, 054012 (2001) [arXiv:hep-ph/0008046]; G. Isidori and 
A. Retico, JHEP 0111, 001 (2001) [arXiv:hep-ph/0110121]; C. Bobeth, T. Ewerth, F. Kruger 
and J. Urban, Phys. Rev. D 64, 074014 (2001) [arXiv:hep-ph/0104284]; S. Baek, P. Ko and 
W. Y. Song, Phys. Rev. Lett. 89, 271801 (2002) [arXiv:hep-ph/0205259]; JHEP 0303, 054 (2003) 
[arXiv:hep-ph/0208112]; T. Blazek, S. F. King and J. K. Parry, Phys. Lett. B 589 (2004) 39 
[arXiv:hep-ph/0308068]; S. Baek, Phys. Lett. B 595, 461 (2004) [arXiv:hep-ph/0406007]. 

[6] S. Baek, J. H. Jeon, and C. S. Kim, Phys. Lett. B 641, 183 (2006) [arXiv:hep-ph/0607113]. 

[7] J. L. Lopez, D. V. Nanopoulos and K. j. Yuan, Nucl. Phys. B 399, 654 (1993) [arXiv:hep- 
th/9203025]; J. L. Lopez and D. V. Nanopoulos, Phys. Rev. D 55, 397 (1997) [arXiv:hep- 
ph/9605359]. 

[8] A. Abulencia et al. [CDF Collaboration], Phys. Rev. Lett. 97, 242003 (2006) [arXiv:hep- 
ex/0609040]. 

[9] S. Aoki et al. [JLQCD Collaboration], Phys. Rev. Lett. 91, 212001 (2003) [arXiv:hep-ph/0307039]; 

A. Gray et al. [HPQCD Collaboration], Phys. Rev. Lett. 95, 212001 (2005) [arXiv:hep- 

lat/0507015]; M. Okamoto, PoS LAT2005, 013 (2006) [arXiv:hep-lat/0510113]; P. Ball and 

R. Fleischer, Eur. Phys. J. C 48, 413 (2006) [arXiv:hep-ph/0604249]. 
[10] Z. Ligeti, M. Papucci and G. Perez, Phys. Rev. Lett. 97, 101801 (2006) [arXiv:hep-ph/0604112]. 

P. Ball and R. Fleischer, Eur. Phys. J. C 48, 413 (2006) [arXiv:hep-ph/0604249]; Y. Grossman, 

Y. Nir and G. Raz, Phys. Rev. Lett. 97, 151801 (2006) [arXiv:hep-ph/0605028]; A. Datta, Phys. 

Rev. D 74, 014022 (2006) [arXiv:hep-ph/0605039]. 
[11] K. Cheung, C. W. Chiang, N. G. Deshpande and J. Jiang, Phys. Lett. B 652, 285 (2007) 

[arXiv:hep-ph/0604223]. X. G. He and G. Valencia, Phys. Rev. D 74, 013011 (2006) [arXiv:hep- 

ph/0605202]. 

[12] M. Ciuchini and L. Silvestrini, Phys. Rev. Lett. 97, 021803 (2006) [arXiv:hep-ph/0603114]; 
M. Endo and S. Mishima, Phys. Lett. B 640, 205 (2006) [arXiv:hep-ph/0603251]; J. Foster, 
K. i. Okumura and L. Roszkowski, Phys. Lett. B 641, 452 (2006) [arXiv:hep-ph/0604121]; 
G. Isidori and P. Paradisi, Phys. Lett. B 639, 499 (2006) [arXiv:hep-ph/0605012]; S. Khalil, Phys. 
Rev. D 74, 035005 (2006) [arXiv:hep-ph/0605021]; S. Baek, JHEP 0609, 077 (2006) [arXiv:hep- 
ph/0605182]; R. Arnowitt, B. Dutta, B. Hu and S. Oh, Phys. Lett. B 641, 305 (2006) [arXiv:hep- 



ph/0606130]; J. K. Parry, Mod. Phys. Lett. A 21 (2006) 2853 [arXiv:hep-ph/0608192]; J. K. Parry 
and H. h. Zhang, arXiv:0710.5443 [hep-ph]. 
[13] A. Lenz, Phys. Rev. D 76, 065006 (2007) [arXiv:0707.1535 [hep-ph]]. 

[14] J. H. Jeon, C. S. Kim, J. Lee, and C. Yu, Phys. Lett. B 636, 270 (2006) [arXiv:hep-ph/0602156]; 
[15] Baek and London in [2]. 

[16] T. G. Rizzo, Phys. Rev. D 59, 015020 (1999) [arXiv:hep-ph/9806397]; K. S. Babu, C. F. Kolda and 
J. March-Russell, Phys. Rev. D 54, 4635 (1996) [arXiv:hep-ph/9603212]. K. S. Babu, C. F. Kolda 
and J. March-Russell, Phys. Rev. D 57, 6788 (1998) [arXiv:hep-ph/9710441]. 

[17] D. Becirevic et al, Nucl. Phys. B 634, 105 (2002) [arXiv:hep-ph/0112303]. 

[18] A. J. Buras, M. Jamin and P. H. Weisz, Nucl. Phys. B 347, 491 (1990). 

[19] J. Charles et al. [CKMfltter Group], Eur. Phys. J. C 41, 1 (2005) [arXiv:hep-ph/0406184]. 

[20] B. Abbott et al. [DO Collaboration], FERMILAB-CONF-97-356-E Presented at 18th International 
Symposium on Lepton and Photon Interactions (LP 97), Hamburg, Germany, 28 Jul - 1 Aug 
1997, and Presented at International Europhysics Conference on High-Energy Physics (HEP 97), 
Jerusalem, Israel, 19-26 Aug 1997 

[21] See Buras, et. al. in [2]. 

[22] Heavy Flavor Averaging Group (HFAG), arXiv:hep-ex/0603003. 

[23] [Particle Data Group], J. Phys. G 33 (2006) 1; B. Aubert et al. [BABAR Collaboration], Phys. 
Rev. Lett. 97, 171805 (2006); K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 98, 181804 
(2007) [arXiv:hep-ex/0608049]; CLEO Collaboration, A. Bornheim et al, Phys. Rev. D 68, 
052002 (2003); B. Aubert et al. [BABAR Collaboration], arXiv:hep-ex/0607106; B. Aubert et 
al. [BABAR Collaboration], Phys. Rev. D 75, 012008 (2007) [arXiv:hep-ex/0608003]; B. Aubert 
et al. [BABAR Collaboration], arXiv:hep-ex/0607096; S. Chen et al. [CLEO Collaboration], Phys. 
Rev. Lett. 85, 525 (2000); K. Abe et al. [Belle Collaboration], Phys. Rev. D 76, 091103 (2007) 
[arXiv:hep-ex/0609006]. 

[24] M. Gronau, D. Pirjol and T. M. Yan, Phys. Rev. D 60, 034021 (1999) [Erratum-ibid. D 69, 

119901 (2004)] [arXiv:hep-ph/9810482]. 
[25] M. Beneke and M. Neubert, Nucl. Phys. B 675 (2003) 333. 

[26] H. n. Li, S. Mishima and A. I. Sanda, Phys. Rev. D 72, 114005 (2005). We thank the authors for 

sending us the updated results. 
[27] C. W. Bauer, I. Z. Rothstein and I. W. Stewart, Phys. Rev. D 74, 034010 (2006). 
[28] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996) [arXiv:hep- 

ph/9512380]. 



